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We simulate quenched QCD with the overlap Dirac operator. We work with the Wilson gauge action at /3 = 6 
on an 18 3 x 64 lattice. We calculate quark propagators for a single source point and quark mass ranging from 
am q — 0.03 to 0.75. We present here preliminary results based on the propagators for 60 gauge field configurations. 



1. Introduction 

The closely related domain wall [1,2] and over- 
lap [3-5] formulations of lattice fermions, with 
their ability to preserve chiral symmetry even at 
finite lattice spacing [6,7], offer an almost ideal 
tool for lattice QCD calculations. Still overlap or 
domain wall fermions can only be implemented 
at a high computational cost. Thus it is im- 
portant to perform exploratory calculations on 
large lattices in order to validate the applicabil- 
ity of these novel formulations. At the same time, 
large scale QCD simulations with overlap or do- 
main wall fermions, because of the benefits of chi- 
ral symmetry, can produce valuable cross-checks 
of observables calculated with more traditional 
quark discretizations, or even permit the evalu- 
ation of observables otherwise out of the reach 
of practical calculations. In earlier studies we 
simulated quenched QCD with overlap fermions 
on a 16 3 x 32 lattice, obtaining results for the 
pseudoscalar spectrum, strange quark mass and 
quark condensate [8,9], for the kaon B parame- 
ter [10,11] and for non-perturbative renormaliza- 
tion constants [9-11]. Although the results of this 
earlier work turned out to be quite satisfactory, 
the size of the lattice, especially its extent in time, 
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were found insufficient for the calculation of some 
observables. For example, the plateaus needed to 
isolate the K mesons in the calculation of Bk ex- 
tended for only two lattice spacings and we could 
not reliably evaluate vector meson masses and de- 
cay constants or baryon masses. This prompted 
us to extend our overlap calculations to a larger 
system. It is important to observe that increasing 
the size of the lattice in simulations with overlap 
fermions does not only entail an augmentation of 
the (already quite heavy) computational cost. In- 
deed, since the calculation of the overlap operator 
requires the use of some suitable approximation 
to the inverse square root of a very large sparse 
matrix, the increase of the size of the matrix may 
produce serious problems of convergence. For 
these reasons, and also in consideration of the 
computational resources available to us, we de- 
cided to use a lattice of size 18 3 x 64, with double 
the extent in time, but only a moderate increase 
in the spatial extent with respect to our former 
calculation. Thus we generated and archived 100 
18 3 x 64 pure gauge field configurations with the 
Wilson gauge action at (3 = 6. We used a 6-hit 
Metropolis algorithm, tuning the acceptance to 
f=a 0.5, performed 11,000 initial equilibrating it- 
erations, after which 10,000 iterations were done 
between each pair of subsequent configurations. 
From a measurement of Wilson loops we deter- 
mined ro/a = 5.36 ± 0.11 for the Sommer scale 
defined by r^F(rQ) — 1.65. With r$ — 0.5fm we 
get a- 1 = 2.11 ± 0.04GeV. 

For each of the above configurations we cal- 
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culated quark propagators with the overlap 
Dirac operator, with p = 1.4, for a single 
point source and all 12 color-spin combinations, 
for am q = 0.03,0.04,0.06,0.08.0.1,0.25,0.5,0.75. 
For brevity we do not reproduce here the for- 
mula for the overlap Dirac operator: we refer 
to [8,9] for the relevant equations and selected 
citations. The results we present here are based 
on the analysis of 60 configurations. For the cal- 
culation of the overlap operator (more properly, 
of its action on a given vector) for the first 55 con- 
figurations we used the Zolotarev approximation 
with 12 poles, after Ritz projection of the lowest 
12 eigenvectors of H 2 . We found, however, the 
Chebyshev approximation to be computationally 
less costly (by approximately 20%) and therefore 
we are now using this approximation. With both 
approximations we impose \D' Dip — x\ 2 < 10 -7 
as convergence criterion. We project out the low- 
est 12 eigenvectors of H 2 also with the Cheby- 
shev approximation and use the evaluation of the 
next highest eigenvalue to set the required de- 
gree of the expansion, which typically varies be- 
tween 100 and 500. Our code has been written 
in F90 with OMP directives and we have been 
running in shared memory mode on 16 and 32 
processor IBM-P690 nodes at Boston University 
and NCSA. 

To conclude this section, we would like to men- 
tion that other large scale quenched QCD calcu- 
lations with overlap fermions have been presented 
in [12-15]. 

2. Light hadron spectrum and quark con- 
densate 

Our first concern has been to verify that the re- 
sults with the new, larger lattice are compatible 
with the values found for the observables in our 
former calculations [8,9] and that the measure- 
ment of the observables can be extended to the 
much larger temporal extent with reliable statis- 
tics. We present here a sample of these checks. 

In Fig. 1 we report the new results for the ratio 

P (t) = G VoAoP (t)/Gpp(t) (1) 

where Gu(t), I J = VqAqP, PP, are non-singlet 
zero-momentum, 2-point functions made from the 
bilinears I and J. On account of the axial Ward 



identity, p(t) should be constant in time. The 
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Figure 1. Results for the axial Ward identity. 

figure shows that p(t) can be reliably measured 
over the entire extent of the lattice for all quark 
mass values and is indeed constant in time. The 




Figure 2. The ratio p of Eq. 1 as function of quark 
mass. 



fit p — A + a(mi + rri2)/ZA + C{am\ + am^) 2 (see 
Fig. 2) gives A = -(0.96 ± 3.35) x 10~ 5 ,Za = 
1.554±0.001 [1.55±0.04],C= (6.87±0.18)xl0- 2 . 
The numbers in square brackets, here and in the 
following, indicate results we found in [8,9]. 

Figure 3 illustrates the value of the effective 
mass from the (0-momentum component of the) 
correlation functions of the pseudoscalar density, 
of the pseudoscalar minus scalar densities and of 
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Figure 3. Effective pseudo-scalar mass plateaus 
for am q — 0.06. 



the temporal component of the left-handed cur- 
rents, for am q = 0.06. We see that the sig- 
nal for the correlation functions extends through 
the whole lattice. These three correlation func- 




Figure 4. Pseudoscalar mass squared as function 
of the total quark mass. 

tions, symmetrized with respect to t — > 64a — t, 
have been fit in the range 12 < t/a < 32, and 
Fig. 4 displays the values thus found for the pseu- 
doscalar mass as well as the linear fit (aMp) 2 = 
A + Ba{pi\ + m-2)j1 which gives 
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PP : A = (6.69 ± 1.76) x 10~ 3 , 
SS : A = (7.27±2.44) x 10~ 3 , 
A = (5.31 ±2.58) x 10~ 3 , 



B = 1.38 ±0.02 
B = 1.37 ±0.02 
B = 1.38 ±0.04 
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Figure 5. Determination of the chiral condensate. 

Figure 5 shows the results for the scalar con- 
densate from the Gell-Mann, Oakes and Renner 
relation, for degenerate quark masses, together 
with the linear fit x — A + Bam q , which gives 
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(1.66 ±0.18) x 10~ 3 
(1.03 ±0.15) x 10~ 2 
(1.67 ±0.24) x 10" 3 
(1.04 ±0.21) x 10~ 2 
(1.17 ±0.27) x 10~ 3 ] 




[PP : B = 1.39 ± 0.03; PP - SS : B = 1.43 ± 0.07] 



Figure 6. Determination of a. 

Finally, we reproduce in Fig. 6 the determi- 
nation of a with the method of lattice physical 
planes: the red line corresponds to the physi- 
cal ratio F K /M K with M K = 0.495 GeV F K = 
0.16 GeV. We find a" 1 = 2.238 ± 0.083 GeV from 
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the PP channel and a" 1 = 2.236 ± 0.098[2.29 ± 
0.09] GeV from the PP-SS channel. 

Figures 7, 8 illustrate new results, which were 
out of reach with a smaller lattice. In Fig. 7 we 
show the quantity log\GBB(t)\ with B = B± = 
(1 ± 7o)e a 6c [v-a C7s<4] u c (a sum over the spatial 
coordinates of B over the time slice at t is im- 
plicit). Since the correlation function with the 
l+7o projector describes forward propagating 
l/2 + baryons and backward propagating 1/2 - 
baryons, while the two parities are interchanged 
with the 1 — 7o projector, in Fig. 7 we plot the 
1 — 7o correlation function against 64a — t. The 
results for the two correlation functions, after this 
inversion of the time axis, are quite consistent, as 
one would expect, and exhibit an extended linear 
region over which one can extract an estimate of 
the lowest baryon masses. The corresponding lin- 
ear fits are also shown in Fig. 7. 



"1+G4,m=0.04" h 
"l-G^rr^O.O*' h 
-0.753'x-11.25 
-O.961"(64-x)-10.72 



1.5 



1 



i 



1 - 




0.5 



I N(940) 

N(1535) 
i/=1/2* 

► /=1/2" 



0.2 



0.4 0.6 

(2m 1 4-m 2 )a 



Figure 8. Baryon masses as function of the total 
quark mass. 



particular the masses in Fig. 8 have been obtained 
with point sources and sinks. While having cal- 
culated point source quark propagators forces us 
to use point sources, we should be able to im- 
prove the statistical significance of the results by 
using suitable extended sink operators. Calcula- 
tions along these lines are in progress. 
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Figure 7. Logarithm of the baryon-baryon corre- 
lation functions. 

In Fig. 8 we collect all of our results for the 
lowest 1/2+ and l/2~ baryon masses. We cal- 
culated the masses corresponding to all combina- 
tions where two of the quarks have the same mass. 
The masses of the l/2 + states appear to fall on a 
single curve, indicating a dependence only on the 
total quark mass. The masses of the l/2~ seem to 
indicate a dependence on the quark mass differ- 
ence also, but the statistical errors are high and 
a more elaborated analysis should be done before 
deriving any such implication from the data. We 
emphasize that these results are preliminary. In 



3. Experimenting with the charm quark 

In the present section we investigate the behav- 
ior of overlap fermions with heavy quarks Q for 
which arriQ <jt. 1. In particular, we are interested 
in understanding whether the charm quark can 
be approximately simulated at lattice spacings of 
around 2 GeV, where the charm has a bare mass 
am c ~ 0.75. Related results, obtained with a dif- 
ferent gauge action, can be found in [16]. 

Though a rigorous investigation of discretiza- 
tion errors would require simulations at a num- 
ber of values of the lattice spacing, at the present 
time we only have a simulation for which a -1 ~ 
2 GeV. We have thus devised a number of tests 
which provide us with some measure of the mass- 
dependent discretization errors present in our 
simulation. They concern a variety of quantities, 
some of which are of direct phenomenological in- 
terest such as the leptonic decays constant of the 
D s meson, f Ds . 
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3.1. Axial and vector Ward identities 

We turn once again to the non-singlet axial 
(AWI) and vector (VWI) Ward identities. In- 
stead of working with the ratio p(t) defined in 
Eq. (1), we consider the quantity 



Z A (t;mi,m 2 ) 



(ami + am 2 ) 
p(t;mi,m 2 ) 



(2) 



with a similar definition for Zy(t;mi,m 2 ) in 
terms of vector currents and scalar densities. 
These quantities are also constant in time, by 
virtue of the AWI and VWI, and reduce to the 
renormalization constants Za and Zy of the axial 
and vector currents in the chiral limit. Za (t) and 
Zy(t) have very nice plateaus in t for all quark- 
mass combinations, though error bars and fluctu- 
ations on Zy(t) are substantially larger. We thus 
fit them to a constant in the time range 14 < at < 
50, thereby obtaining mass-dependent "renormal- 
ization constants" ZA(mi,m 2 ) and Zy(mi,m 2 ). 
The fits to Za{£) are all excellent, except for the 
mass combination ami = am 2 = 0.75 where there 
is some evidence that the signal is lost for a hand- 
ful of points around the center of the lattice (cf. 
Fig. 1). The situation is much less satisfactory 
for Zy (t) , where the fluctuations make it difficult 
to obtain reliable fits. 

The quantities Z^y^mi, m 2 ) are of particular 
interest here, because their mass dependence is 
solely due to discretization errors. Given our defi- 
nition, the leading mass dependence is of the form 
a 2 (mi + 777,2) and a 2 (mi — 7712) 2 / (mi + 7712). In 
Fig. 9 we plot ^(mi, m 2 )/^(0, 0) as a function 
of (ami + 0777,2) for all mass combinations, were 
Za(0, 0) is obtained as described below. We con- 
sider Zy (777,1,7712) only for (ami + 07772) < 0.2, 
since it is not clear that we can reliably fit Zy (t) 
for larger masses. The values of ZA(mi,m 2 ) for 
(ami + 07712) < 0.2 display linear behavior, indi- 
cating that mass-dependent discretization errors 
are dominated by a 2 (mi +m2)-type terms in this 
range of masses. Therefore, in this mass range we 
fit ZA{m q , m q ) to: 

Z A (m q ,m q ) = Z A (0,0)(l + 2(aB)(am q )) , (3) 

where we only consider degenerate combinations 
to eliminate possible a 2 (mi — m 2 ) 2 /(mi + m.2) 



terms. We find Z A (0,0) = 1.5544(15), which is 
entirely compatible with the value given after Eq. 
(1), and aB = -0.102(5). 
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Figure 9. ZA,v(m 1 ,m 2 )/ZA(0,0) vs. (ami + 
07772). The straight line corresponds to the fit 
described around Eq. (3). The two filled symbols 
approximately correspond, from left to right, to 
a D s and an rj c meson. 

The first point worth noting is that the coef- 
ficient B of the a 2 (mi + 7772) discretization error 
is small, B ~ 225 MeV. The second interesting 
feature is that deviation from linearity becomes 
statistically significant (i.e. larger than 3 standard 
deviations) only for (ami +am 2 ) > 0.28. We are 
actually able to consistently fit .£4(777!, 7772) all 
the way up to ami = 07772 =0.5 with the addition 
of two terms proportional to a 3 (mi + m2) 2 and 
a 2 (mi — m 2 ) 2 / (mi + m 2 ), and all the way up to 
ami = 07772 = 0.75 with an additional term pro- 
portional to a 4 (mi +77i2) 3 , though the coefficients 
obtained in the latter fit are only marginally con- 
sistent with those from our two other fits. In these 
fits to extended mass ranges, we further find that 
the coefficients of the higher order terms are ap- 
proximately one half to a full order of magnitude 
smaller than aB. The net result is that discretiza- 
tion errors inferred from ^(roi, 7772) are approxi- 
mately 10% at (ami,am 2 ) = (0.06, 0.75) and 15% 
at (ami, 07772) = (0.75,0.75), mass combinations 
which corresponds roughly to a D s and an rj c 
meson, respectively. Moreover, chiral symmetry 
seems to be good, as evidenced by the compat- 
ibility of ZA(mi,m 2 ) and Zy(mi,m 2 ) at small 
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masses. Thus, there is no evidence for a major 
breakdown of overlap fermions up to am, = 0.75 
and (ami + 0,1712) = 1-5. 

3.2. 2-point functions, pseudoscalar meson 
masses and decay constants 

We have performed fits to the zero-momentum 
PP, PP-SS, A P, PA a , A a A and J L J L corre- 
lators, all of which have the pseudoscalar meson 
as the lightest contributing state. The fit form 
used is: 



£<J(aOJ(0)) 

X 

2M U 



e -Muxo _|_ Slje -Mu(T-xo) 



(4) 



where the fit parameters are Zu and M/j, and 
sij = — 1 for IJ = AP,PA and +1 otherwise. 
For all of our quark masses, the fits are excel- 
lent. For each correlation function, we define the 
corresponding pseudoscalar decay constant: 




AP,PA 



where, in our conventions, f v — 131 MeV. The 
point of considering these three definitions is that 
they differ in their discretization errors. In fact, 
we expect these differences to be at most equal 
to the mass-dependent discretization errors iden- 
tified in ^(mi, m 2 ). Indeed, fp P /fp A is simply 
ZA(mi,m2)/ZA.(0,0) to the extent that the fits 
to the continuum-limit parametrizations of Eq. 
(4) do not modify the relative discretization er- 
rors in the amplitudes Zu; that Vo applied on 
GA p(t) yields a value of the pseudoscalar mass 
which is the same as the one determined by the 
fits; that finite-volume effects do not distort the 
various correlation functions differently. And un- 
der these same conditions, all other ratios of de- 
cay constants should be less than or equal to 
fp P /fp A - Moreover, the correlation functions 
are affected by topological zero modes differently. 
In particular, PP— SS and JqJq are free of zero 



modes. Differences at small quark masses thus 
give a measure of zero-mode contamination. 

In Fig. 10 we plot fp J versus M/j, for all 
IJ. No statistically significant difference in 
the masses obtained from the various correlation 
functions is observed. At low quark mass this in- 
dicates that there is no evidence for zero-mode 
contamination in our fits. At high quark mass, 
the agreement is reassuring but not significant as 
discretization errors on the mass are determined 
by the action, which is the same for all corre- 
lation functions. That is not the case for the 
decay constant, where different correlation func- 
tions will induce different discretization errors. 
Indeed, we find that the difference increases as 
the quark mass is increased. For the U D S " , we 
find a variation of approximately 10% and it is 
15% for the "?7 C " , which is entirely consistent with 
the results obtained from Za (mi, 7712), as antici- 
pated. While significant, these variations do not 
indicate a major breakdown of the theory, even at 
ami+ani2 = 1.5. This is further confirmed by the 
fact that with a charm quark such that the corre- 
sponding "£) a " meson has a mass <~ 1.95 GeV (i.e. 
20 MeV below experiment), our "r/ c " has a mass 
- 2.93 GeV which is only 50 MeV (i.e. 2%) below 
experiment. Moreover, the preliminary value of 
f Ds that we obtain is 267(15) (24) MeV, which is 
in good agreement with the most recent quenched 
calculation [17] and experimental measurement, 
f Ds = 285(19) (40) MeV [18]. 

3.3. Charm quark summary 

We have simulated the charm quark with over- 
lap fermions on a lattice with a" 1 ~ 2 GeV. We 
have investigated the effects of mass-dependent 
discretization errors on the axial Ward identity as 
well as on the pseudoscalar decay constant. We 
find that these effects are relatively small: 10% for 
the "L> s " and 15% for the V'. 3 We further find 
a value of fo 3 which is in good agreement with the 
latest quenched and experimental results. Thus, 
there is no sign of a major breakdown of the for- 
malism even for bare quarks masses am q — 0.75. 
Of course, these preliminary findings require fur- 



3 It should be noted, however, that these figures come from 
ratios of quantities in which some discretization errors may 
cancel. 
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Figure 10. fp vs. Mp as obtained from vari- 
ous 2-point functions. For clarity, only results for 
mesons composed of degenerate quarks are plot- 
ted, with the exception of our fiducial D s meson. 
The vertical lines correspond, from left to right, 
to P = D s and P = r) c . 

ther investigation and, in particular, the apparent 
smallncss of discretization effects should be con- 
firmed by simulations on finer lattices. 

4. Weak matrix elements 
4.1. Bare matrix elements 

We have calculated the K°-K° mixing parame- 
ter Bk following the procedure detailed in [10,11]. 
B>k is defined as 



(K \O AS = 2 (ri\K°) 



with ( 7 ^ = 7m (1 T75 )) 



in 

OaS=2 



(6) 



(7) 



The bare matrix element (K°\Oas=2\K°) is ex- 
tracted from the 3-point function 

Cjoj(xo,yo) = Y,( J o^)O b AsU(0)Jo(y)) (8) 



and the 2-point function 



(9) 



which are both free of unphysical zero modes [19]. 

The large temporal extent of our lattice al- 
lows us to choose ranges 10 < xq/cl < 20 and 



44 < yo/a < 54 for the kaon source and sink, 
which are both safely in the asymptotic region 
Xq ^> a, yo <C T and for which backward contri- 
butions due to the toroidal lattice geometry are 
negligible. As a consistency check, we have ver- 
ified that within these ranges the pseudoscalar 
meson mass extracted from Cjoj(xo,yo) agrees 
very well with the ones extracted from various 
2-point functions (cf. Section 2). 

We have also computed the bag parameters 

3/2 

By s of the electroweak penguin operators in the 
AI = 3/2 channel of the AS* = 1 effective Hamil- 
tonian. Following the convention of [20], we de- 
fine 



B^ 2 (n) = lirn^ 2 



(7r+\0 3 7 /2 (v)\K+) 



m^o |( 7 r+|P(/i)|0)(0|P( )U )|^ 
{^ + \0T^)\K + ) 



B%'{n) = lim 



m^o 2{ir+\P{p)\Q){Q\P{^)\K+) 



(10) 



(11) 



3 /2 

where 7 is given by 

3 7 /2 =srfd(u 1 K U -drfd) + srfuurfd (12) 

3/2 

and Og is the corresponding color-mixed oper- 
ator. We work in the SU (3)y limit m u = = 
m s , in which all eye contractions cancel. 

In contrast to B^ , finding appropriate sources 
to cancel all zero-mode contributions is not 

3/2 

straightforward for By s due to the L-R struc- 
ture of the operator. For this preliminary study, 
we therefore choose to simply report results ob- 
tained with pseudoscalar sources on the 3-point 
function 

C 7/8 (x ,y ) =^(P(x)(O? / / 8 2 ) bare (0)P(y)> (13) 



and the 2-point function 



Cpp(xo) = £<P(a:)P(0)> 



(14) 



and postpone a more careful investigation of this 
issue to a future publication. 4 

We identify nice plateaus in Cy$(xo, yo) in the 
range 18 < xo/a < 23 and 41 < yo/a < 46. In 

4 Pseudoscalar sources are not the optimal choice, since 
they allow for 1/m 2 divergences in the 3-point functions. 
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this region, the mass extracted from the 3-point 
function agrees very well with the mass extracted 
from various 2-point functions (cf. Section 2). 
4.2. Non-perturbative renormalization 

We perform all renormalizations non- 
perturbatively in the RI/MOM scheme following 
[21]. Thus, we fix gluon configurations to Lan- 
dau gauge and numerically compute appropriate, 
amputated forward quark Green functions with 
legs of momenta p. The renormalization of 
is detailed in [11]. In order to renormalize By^, 
we calculate the renormalization constants of 
a full set of dimension 6, 4-flavor operators in 
the basis suggested by [22]. In this basis, O^ 2 

renormalizes as and O^ 2 as — 2Qf. These 
2 operators mix among each other but chiral 
symmetry protects them from mixing with other 
operators. The renormalization condition vio- 
lates chiral symmetry at low p 2 , but we found, 
that the violation is tiny. 

The renormalization pattern for the B param- 
eters is given by 



5 10 15 5 10 15 20 



z 7 M 






Figure 11. The coefficient matrix IZf* vs. fi. One 
can see, that the data are very well described by 
Eq. (17). 



3/2 
7 

3/2 




(15) 



In order to determine the Z^(fi) we define 



TZgfy) = lim 



Tr(P s r s (p 2 ,m)) 2 



Tr^r^m)) 



(16) 



where Ti(p 2 , m) is the amputated Green function 

3/2 

of the operator O i and Pi the corresponding 
projector (cf. [22]). After performing an appro- 
priate chiral extrapolation, we must isolate the 
"perturbative part" of this ratio to get the renor- 
malization constant. In order to do this, we fit 
the matrix 7Z m (/i) to the form 



An B 
-T + - 

x {an) 2 - 

RI/ 



(/*) 



(17) 



where the running of Z^j (n) is given by NLO 
continuum perturbation theory. The l//i 2 and 
l/// 4 terms are motivated by the OPE in 1/p 2 



of Eq. (16), while the term cx (a//) 2 takes into 
account leading order discretization effects. 

In Fig. 11 the 1lff{Li) are shown together 
with fits of the form Eq. (17) in a range fi G 
[y/2, 3] GeV. The fit describes the data very nicely 
even far outside the fit range and is stable against 
addition of extra terms and variations in the fit 
range. Some of the fit parameters have numeri- 
cally found to be consistent with and have been 
eliminated. 

4.3. Physical results 

In Fig. 12 we show the renormalized value of 
B K in the MS-NDR scheme at 2 GeV. We have 
performed both a linear and a chiral log fit of the 
form 



B K = B 1 -6 



(M_\\M 
UttFJ A b 



(18) 



At the physical point, the result of both fits 
is virtually indistinguishable and we obtain 
B MS-NDR(2GeV) = 0.61(7), where the error is 
statistical only. This result is compatible both 
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Figure 12. B% s " NDR (2GeV) vs. the pseudoscalar 
meson mass. The continuous line gives the result 
of a linear fit and the dotted line of a chiral log 
fit of the form Eq. (18). 



with our earlier result from smaller lattices [11] 
and with a recent determination by the MILC 
collaboration [15]. 

In Fig. 12 we plot the renormalized value of 
B^l in the MS-NDR scheme at 2 GeV as a func- 
tion of Mp. We postpone the study of this mass 
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Figure 13. (s^f) 
doscalar meson mass 
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(2 GeV) vs. the pseu- 
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